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We study the classical and quantum dynamics of generally co- 
variant theories with vanishing Hamiltonian and with a finite 
' number of degrees of freedom. In particular, the geometric 

meaning of the full solution of the relational evolution of the 
£NJ , degrees of freedom is displayed, which means the determina- 

tion of the total number of evolving constants of motion re- 
I quired. Also a method to find evolving constants is proposed. 

The generalized Heinsenberg picture needs M time variables, 
■ as opposed to the Heisenberg picture of standard quantum 

mechanics where one time variable t is enough. As an appli- 
cation, we study the parameterized harmonic oscillator and 
the SL(2, R) model with one physical degree of freedom that 
j 1 ^ ' mimics the constraint structure of general relativity where a 

Schrodinger equation emerges in its quantum dynamics. 



x. 

Key words: Evolving constants of motion. 



INTRODUCTION 



Describing evolution of the degrees of freedom of generally covariant theories is an 
unsolved puzzle, and constitutes one of the challenges of the human thinking of our 
time. The study of generally covariant theories has been motivated by general relativity, 
which has this peculiar property (see for instance fl|]). In gravity, general covariance 
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means the theory is diffeomorphism invariant, and this symmetry of gravity implies at 
the Hamiltonian level that the theory has not a genuine Hamiltonian for describing the 
evolution of the degrees of freedom of the gravitational field, rather, dynamics is gauge; 
generated by the first class constraints of the theory. This is the so-called problem of 
time in classical general relativity |2j . 

On the other hand, if the classical regime of general relativity is only a limit, which 
emerges in a suitable way from its fundamental quantum behavior , then the theory is in 
trouble. Our standard methods of quantization crash and do not apply to the particu- 
lar physical situation raised by general relativity. Standard quantization methods in field 
theory are background-dependent, quantum gravity needs a background independent pro- 
cedure in its quantization. So, how to make compatible the symmetry of general relativity 
with the quantum theory? what is the meaning of evolution in quantum gravity? Loop 
quantum gravity answers the first question, because the quantization of the gravitational 
field is carried out in a background independent way ||. With respect to the second, it 
remains as an open question. Among the several proposals for describing the evolution 
of the degrees of freedom of generally covariant theories, I find Rovelli's proposal as one 
of the most creative ones J§-§|. 

Here, following the spirit of relationism, which is the heart in Rovelli's point of view, we 
analyse the 'problem of time' in generally covariant theories with vanishing Hamiltonian 
and with a finite number D of degrees of freedom. To obtain the relations involving the 
coordinates and momenta of the unreduced phase space 7 ez with the physical states that 
label the points of the physical phase space 7^, we need to start from the embedding 
equations which give the dependence of the coordinates and momenta with respect to 
the M time variables t m as well as the 2D physical states {q a ,p a )- These equations 
constitute the classical version of the generalized Heisenberg picture, which arises when 
these equations are promoted to quantum operators in the reduced Hilbert space of the 
theory. By plugging the expressions of the time variables t m in terms the original canonical 
variables into the expressions of coordinates and momenta, we get the full relational 
evolution of the phase space degrees of freedom for any physical state of 7^. This way 
of expressing the full solution of the dynamics of generally covariant theories constitutes 
the full set of evolving constants of motion required in their dynamics, and is displayed in 
Sect. I. In addition, an alternative mechanism which generates also the evolving constants 
of motion is proposed. Of course, we study also the quantum version of the evolving 
constants of motion. In sect. II, we analyze the parameterized harmonic oscillator (as an 
example of parameterized systems). In Sect. Ill, we continue the study of the SL(2, R) 
model which constraint algebra mimics the algebra structure of general relativity. Due 
to the fact a Schrodinger equation emerges in its quantum dynamics, we compare the 
generalized Heisenberg picture (related with the evolving constants of motion) which 
needs M time variables with the Schrodinger picture which singles out one time variable 
only. We also especulate on the classical limit generally covariant theories and its possible 
relation with the full set of evolving constants of motion. Our conclusions are summarized 
in Sect. IV. 
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I. RELATIONAL CLASSICAL AND QUANTUM DYNAMICS 



Classical dynamics. The classical dynamics of a constrained theory with a finite 
number of degrees of freedom characterized by first class constraints is as follows H . The 
theory is obtained from the Hamiltonian form of the action 

S[q\Pi, X m ] = jH dr \^- Vi - X m C m {q\ Pi )\ , (1) 



which is invariant under arbitrary reparameterizations of the parameter r. Hence, r is 
a non physical coordinate time. The unreduced phase space 7 ea; is coordinated by the 
canonical pairs (q\Pi) \i = 1,2,..., N. The canonical 2-form on j ex is u ex = dpi A dq\ 
Thus, (7,o>) is a symplectic space. The variation of the action S[q l ,pi,X m ] with respect 
to the canonical coordinates q\pi gives the equations of motion 

dtf = } m dC m{.q\Pi) 

dr dpi 

dpi xm dC m {q\pi) 



while the variation of the action with respect to the Lagrange multipliers X m gives the 
constraint equations 

C m = C m (q i ,p i ) = 0, m = l,2..,M. (3) 



The variation of the action has been done under the standard boundary conditions 
q l {r s ) = ql; s = 1,2, namely, the allowed paths are those with fixed values for the con- 
figuration variables at the boundary points r s . The boundary conditions can be changed 
and thus to modify the action by suitable boundary terms to allow the gauge symme- 
try generated by the constraints ||. The constraints generate Hamiltonian vector fields 
Xdc m , which are tangent vectors to the constraint surface, given by 

X dC - dCm - + dCm - (4) 
dpi ()(/' i)(f dpi 



More important, the integral curves of these Hamiltonian vectors fields constitute the 
gauge submanifold or the orbits of the constraint surface, and the dynamics of the system 
with respect to r is the unfolding of this gauge symmetry, i.e., dynamics is gauge. 

The first class constraints satisfy, in general, a non Lie algebra 

{C m ,C n } = C mn l (q t ,Pi)Ci , (5) 
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and the number of independent physical degrees of freedom of the theory is D = N — M. 
The constraint surface defined by the constraint equations (y) is a (2D + M) -dimensional 
manifold. The constraint surface can be parameterized by the set of independent coordi- 
nates (q a ,p a , t m ), where {q a ,p a ), a = 1,2, D are (local) canonical variables which coor- 
dinatize the open sets of the physical phase space j p h of the theory, and t m ,m = 1,2, M 
coordinatize the orbits, i.e., the gauge submanifold of the constraint surface generated by 
the first class constraints. Notice that the canonical coordinates qo a , and po a are the 
physical observables of the system. Of course, they satisfy {q a ,pb} = o~i on the physi- 
cal phase space, and the symplectic form on ^i ph in these coordinates is u ph = dp a A dq a . 
Therefore, the general solution of the dynamics of the constrained theory can be expressed 
as 

q l = q l (t m ;q a ,Pa), (6) 
Pi=Pi{t m ;q a ,p a ). (7) 

It is important to emphasize that that such dependence is local. For instance, in the 
case in which the physical phase space 7^ is compact, a finite set of physical observables 
(q a ,Pa) is needed to coordinate the open sets of 7^ due to its compactness. So, the 
constraint surface looks like a 'fibre bundle' P(j p h, Orbits), the constraint surface being 
the total space P, the physical phase space "f p h being the base space, and the orbits being 
the fibers of the bundle. In the generic case, P{^ p h, Orbits) is locally trivial. This means 
that non global gauge condition is allowed in general, and that local gauge conditions 
associated with each open set of the physical phase space can be specified only. 

At the same time, the full solution of the theory implies to give the dependence of the 
physical observables q a , and p a of the system in terms of the coordinates of the unreduced 
phase space 

? = ?tf,Pi), (8) 

Pa=Pa(q\Pi) , (9) 

as well as the orbit coordinates t m 

t m = t m (q\ Pi ). (10) 

What these equations tell us is that one single internal time variable is not enough to 
describe the evolution of the system, rather, M internal time variables are needed. In the 
way the full solution has been expressed in (0), and (|^), these M time variables are t m , 
m = 1, 2, M. Notice also that these M time variables are internal clocks, given by (PH|). 
One of the properties of these internal clocks t m is that they do not run taking increasing 
values of t m when time goes on. In fact, they can run 'forward' and 'backward' depending 
on the values of the coordinates and momenta the system is reaching through Eq. (|10|). 
Other property is that these clocks can run with different 'speeds' for the same reason. 
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So, the meaning of time that arise in generally covariant theories is completely different 
with respect to the monotonous function we are familiarized with. In || was showed that 
Eqs. (^j and (|7|) can be obtained by a combination of a canonical transformation plus 
Hamilton- Jacobi techniques. That approach implies the modification of the original set 
of first class constraints. This can be done, but it is not necessary in principle. Moreover, 
the full solution requires (|T0|) (missing in Ref. |§) as we have seen, and more important, 
it is the combination of Eqs. (0), and (|7]) with Eq. ( [TDD which leads to the relational 
description of the dynamics of the system, as we will see it later on. 

It is worth to mention the relationship between the time variables t m and the full 
gauge transformation generated by the first class constraints C m . Assuming that the full 
gauge transformation of the original canonical variables is given by 

v'i=p'i{q\Pi^ m {r)), (11) 

with a m the m gauge parameters involved in the gauge transformation. Then, by plugging 
(|) and (|7|) into the right hand side of (|TT|), we get 

q' l = q\t' m -~q\Pa), (12) 
p\=Mt' m ;q a ,Pa), (13) 

where the functional dependence in the right-hand side of the above equations is exactly 
the same as that given by ©-(0) but with 

t' m = t' m (t m ,a m ), (14) 

which relates the time variables t m after and before of any finite gauge transformation of 
the canonical variables (|TT|). 

The map : P(^y ph , Orbits) -> j ex , (j)(q a ,p a ,t m ) i-> (q l (t m ;q a ,p a ),pi(t m ;q a ,p a )) allows 
us to define on the constraint surface P^ph, Orbits) the pull back Q = (fi*u ex of the 
canonical form u ex = dpi A dq l on ^y ex , which is degenerate. Thus, the geometry of con- 
strained systems involves three spaces: the unconstrained phase space (jex^ex), the con- 
straint surface (P(j p h, Orbits), Q = <p*u ex ), and the physical phase space (y p h,Wph)- The 
map that connects the constraint surface and the physical phase space is the projection 
7r : P(j p h, Orbits) — > j p h, Tr(q a ,p a ,t m ) i— > (q a ,p a )- Due to the fact that the Hamiltonian 
vector fields @ are tangent vectors to the orbits, they can be expressed in terms of the 
local coordinates q a ,p a ,t m of the constraint surface as 

X dCm = {C m ,t n }(q a ,Pa,n^- (15) 

The observables (§), @, and the orbits (|T0|) also generate Hamiltonian vectors fields, 
their restriction on the constraint surface are 
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d „ d 
d d 

x dtm = {t m A a }^ a + it m M^- a + i tm ^}^ > ( 16 ) 

where it is understood that all the quantities are evaluated in the point (q a ,p a , t m ). Thus, 
{X d c m , X d ga, X d p a } is a basis, naturally adapted to the involved geometry, of the tangent 
space of the constraint surface. The vectors X dCm play the role of vertical vectors because 
they have vanishing projection on the tangent space of j p h, d7r(X dCm ) = 0. X d qu , and 
X d p a are the horizontal lifts on the constraint surface of the coordinate basis on 7^, 
G?7r (X d ga) = jj- , drt (X d p a ) = —^=5-. In summary, the solution of the dynamics of the 
constrained system means to specify Eqs. (|B|)-(|rO|). This fact, rises a new problem: 
the problem of the meaning of physical time of generally covariant theories, i.e. the 
specification of an internal clock in the framework of the theory with respect to which to 
describe the evolution of the degrees of freedom of the theory in a gauge invariant way. 
Let us explain, the dynamics with respect to r is given by 

q l (r) = q\t m (r);q a ,Pa), (17) 
p t (r)=p l (t m (r) ] q a ,Pa), (18) 



for any physical state (q a ,p a ) of the system. So, this dynamics is non gauge- invariant, 
i.e., it depends on r. The question is, can we describe evolution of the system in a gauge 
invariant way? The answer is yes. At first sight, this sounds strange or impossible in 
a system with gauge freedom. To see how this can be achieved, let us plug the time 
variables (]T0|) into the full solution @, and (0) 

q l = q l (t m (q\ Pi );q a ,Pa), (19) 
Pi=Pi{t m {q\Pi)-,<t,Pa)- (20) 

Last equations are very important, they relate the original phase space variables q l , and 
Pi with the physical states of the physical phase space (q a ,p a )- These equations admite 
two, related, interpretations. First, they give the relational evolution of the coordinates 
q l and the momenta Pi for any fixed point (q a ,p a ) of the physical phase space, i.e., it 
is possible to choose M coordinates denoted by q m (or momenta p m ; or a combination 
of both) as 'clocks' and describe the evolution of the remaining set of coordinates and 
momenta as functions of the q m for any physical state (q a ,p a ) of the system. Second, if we 
fix the values of this M coordinates, say q m = q* m then, the before mentioned expressions 
of coordinates and momenta give M-parameter families of physical observables defined 
on 7 P h, q* m being the parameters. Eqs. (|T9"D, and (|20"D are evolving constants of motion 
in the sense of Rovelli This concept captures the essence that the before mentioned 
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observables (defined on 7^) describe the relational evolution of the coordinates q % and 
momenta pi, and at the same time they are physical observables. 



Let us consider particular cases of (fTO), say 



q m , m = l,2,...,M, (21) 

then (|T9|), and ( p0|) acquire the form 

q m = q m , m = l,2,...,M, (22) 

q l = q l (q m ;q a ,Pa), i = M + 1, N , (23) 

p l =p l (q m ;q a ,Pa), i = l,...,N. (24) 

Thus, the 'clocks' are given by q m and last two pairs of equations are the evolving constants 
of motion involved. Other particular case is given by 

t m = p m , m = l,2,..,M, (25) 
and (|19D , and (p0|) acquire the form 



q l = q l (Pm;q a ,Pa), i = l,2,...,N, (26) 
Pm=Pm, m = l,2,...,M, (27) 

Pi=Pi(p m ;q a ,p a ), i = M+l,...,N. (28) 

In this case, the 'clocks' are p m and last two pairs of equations are the evolving constants 
of motion required. 

As we have seen, the general relations involving the coordinates q l and momenta Pi 
with the physical states (q a ,p a ) is given by (0), and (EDI). The explicit form of (|T9D , and 
(pQ|) could be complicated for particular theories, but this fact would rise technical rather 
than conceptual difficulties (see |IU| , P^T| for the opposite viewpoint where the authors rise 



questions on interpretation, consistency, and the degree to which the resulting quantum 
theory emerging from the before classical dynamics coincide with, or generalizes, the usual 
non-relativistic theory). Thus, Eqs. (0), and (pOQ constitute the full set of evolving 
constants needed in the relational description of the dynamics of generally covariant 
theories with a finite number D of degrees of freedom. The solution fll9|), and ( |20| ) sits 
in the spirit that in covariant theories there is non privileged observable with respect to 
which to describe evolution, and that only relational evolution makes sense. From this 
point of view, general covariance forces us to use relational evolution, namely, to describe 
the change of some variables of the system with respect to the others. This is the essence 
of relationism, which appears to be the natural language for describing the evolution of 
the degrees of freedom of generally covariant theories @| 



7 



In addition, in this paper, we propose an alternative mechanism to generate the evolv- 
ing constants. This mechanism is essentially to compute the action of the Hamiltonian 
vector fields Xac m on some evolving constant E 1 

X dCm (E) =: E m . (29) 

The evolving function E 1 depends on the canonical coordinates of the unconstrained 
phase space q\ and pi as well as on the canonical coordinates q a , and p a of the physical 
phase space. Therefore, in the computation of the action of the Hamiltonian vector fields 
(f|) on the evolving function we can proceed in two ways. First, taking the observables 
q a , and p a constants in the dependence of the evolving function E. This can be done 
because q a , and p a are constant along the orbits, b) Taking the explicit dependence of 
the physical observables in terms of the canonical variables of the unconstrained phase 
space given by @, and ([5]). Of course, both approaches lead to the same results. The 
repeated application of the Hamiltonian vector fields on the new evolving constants E 2 , 
gives another evolving constants, and so on until no new evolving constants are 
obtained, and the process ends. From the knowledge of the evolving constants and the 
expressions of the physical observables, the full solution of the dynamics of the system 
encoded in Eqs. (§)-(E!3) * s obtained. 

Quantum dynamics. Let us begin with the quantum description of the system. We 
use the Dirac method. In the same way as in the classical dynamics we have three spaces 
(jex, w ex ), (P{y p h, Orbits), fl), and (j p h, 0J p h). In the quantum theory we have three Hilbert 
spaces; the unconstrained Hilbert space 7i or a suitable extension of it if the constraints 
have continuum spectrum, the physical Hilbert space Ti. p h ys , and the reduced Hilbert 
space TC r obtained by projecting H p h ys - Suppose we have solved the quantum theory in a 
full way, i.e., we have the physical Hilbert space TC p h ys of the theory. A general physical 
state | 4>) of the system is killed by all the constraints of the theory C m | ■0) =| 0), and it 
is given by 

\n 1 ,n 2 ,...,n D ). (30) 

ni,U2, —,nr> 

in Dirac notation. Here, the physical states are labelled by the quantum numbers n a , 
a = 1,2, ...,D which come from a complete set of commuting physical observables O a , 
a = 1, 2, ...,D of the system 

O a | ni,n 2 , n D ) = 0(n a ) \ ni,n 2 , n D ) . 

Of course, these quantum observables are combinations of the physical observables (f, 
and p a . We have come to the heart of the problem, how to describe relational evolution 
in the quantum theory. 

Quantum evolving constants. Let us see how the quantum version of the classical 
evolving constants looks. The idea is to search for a representation of the physical states 
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(j30p in the reduced Hilbert space associated with the physical phase space of the system. 
Explicitly 

ij(q a ) = (q a | ifj) = c nun2 _ nD (q a \ n h n 2 , ...,n D ) . (31) 

ni,ri2,...,no 

The inner product in the Hilbert space 

<lH0 = Jdfi(q a )r(Q a mq a ), (32) 

can be determined with the condition that the operators q , and p a be hermitian operators 
and with the implementation of the action of the operators q , p a on this Hilbert space. 
Notice also that is always possible to build creation and annihilation operators a a = 
q + ip a , a\ = q — ip a for each pair of canonical operators q , and p a because the number 
of these operators is even. a a , a\ can help in the construction of TC r . 

With the before machinery, the quantum version of the evolving constants is 

$ = q\t m {q\ Pi )-f3a), (33) 
Vi=Pi{t m {q\Vi)-,T3 a ), (34) 



or, equivalent ly, 

(1> 1 5* I ^ = (1> I q l (t m ( q \ Pl y,T3a) I V»> , (35) 

(i> | Pi | 1>) = (r/> I Pi{t m {q\pi)Cq 3 a ) I V> , (36) 

where the mean values are computed with the inner product . In the case of parame- 
terized systems, last equations reduce to the standard ones which describe the evolution 
of the position and momenta operators as well as the evolution of the mean values of the 
position and momenta operators in the Heisenberg picture. Of course, the well-known 
ordering problems for the operators might appear here too. 



II. PARAMETERIZED HARMONIC OSCILLATOR 



A. Classical dynamics 



In order to make these ideas concrete, let us consider a familiar example: the param- 
eterized harmonic oscillator, which action is 
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S= dr 



dx dt ( p 2 1 2 2 N 

— p + —p t - A [pt + h -ww x 

dr dr \ 2m 2 



(37) 



The unconstrained classical space T is coordinatized by the canonical pairs {x,p), and 
(t,p t ). By doing the variation of the action with respect to x, p, t, and p t we find the 
equations of motion 

^ = -AW*,^ = A^,^ = 0,# = A. (38) 
ar ar m ar ar 

The variation of the action with respect to the Lagrange multiplier A gives the first class 
constraint 

C = p t + ^- + \mu 2 x 2 . (39) 
2m 2 

The classical dynamics is the unfolding of the gauge symmetry of the system generated 
by the first class constraint C. The gauge orbit on the constrained surface C = is the 
integral curve of the Hamiltonian vector field 

X *c = -i--^r + mu 2 x?-. (40) 
ot m ox op 



If we have a solution x(r), p(r), t(r), and Pti^r) of the equations of motion (|38|) , any 
other solution x'(r), p'(r), t'(r), and pt'(r) can be found through the relations 

x'(t) = cos (9(t))x(t) H sin (9(t))p(t) , 

p'(t) = — muj sin (6*(r))a;(r) + cos {6{r))p{r) , 

a; 

P/ = , (41) 

that connect all the solutions, while the Lagrange multiplier transforms as 

A'(r) = A(r) + ^(r) , (42) 
a; 

in order to leave the action invariant , here 9(r) = ^r-^- 

Let us construct the general solution in a given gauge. We choose the gauge A = 1. 
We still have one gauge fixing to impose at r = 0. We choose t(0) = 0. Using the 
constraint equation and the solution of (|38l), we obtain 
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x(t) = A cos (ur) + B sin (ojt) , 
t(r) = r 1 

p(r) = —moo A sin (cur) + muB cos (cjt) , 

Mr) = -\mu\A 2 + B 2 ) , (43) 



where (A, 5) are the physical observables that coordinatize the physical phase space of 
the system, which is R 2 . It is clear that x, t, p are non-observables (they depend on r). 

The two physical observables (A, B) can be expressed in terms of the phase space 
variables as 

A = cos (u>t)x sin (ut)p , 

mui 

B = cos (ujt)p + sin {ut)x (44) 



Notice that A = x(t = 0) = Xq, and B = p ^ ~ ' = i.e., the position x of the 
harmonic oscillator when the internal clock measures t — 0, and the momentum p when 
the internal clock measures t — are (physical) observables. Moreover, Eq. ([H]) means 
that the precise combination of the position x = X and the momentum p = P of the 
harmonic oscillator, when the internal clock indicates t = T in the form expressed by the 
formula ( [4*4] ) is an observable of the (composed) system: harmonic oscillator + internal 
clock. These observables have vanishing Poisson brackets with the first class constraint C 
as required by the formalism of constrained systems. Actually, the Dirac method requires 
observables to have weakly vanishing Poisson brackets with the first class constraints. 
Here, the observables A, B have strong vanishing Poisson brackets with the constraint 
C. The Poisson brackets between A and B in the physical phase space reads 

{A,B} = —. (45) 
mui 



Classical evolving constants. From (^), we obtain the evolving constant 

x = xq cos (ut) + sin (ujt) , (46) 
moo 



of the system. As before mentioned, last equation admits two, related, interpretations. 



First, for any fixed point (A, B) (equivalently (xo,p )) of the physical phase space, ([46]) 
gives the relative evolution of the configuration variables x, and t of the system 

x = X(t; x ,p ) = xq cos (ut) + sin (cut) . (47) 
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Second, for any fixed t, it gives a one-parameter family of physical observables, t being 
the parameter, on the physical phase space. 

Generation of evolving constants. We define the function E l on V 

E 1 (x, t,p,p t ) : = x — Xq cos (ut) — — sin (ut) . (48) 

mu 

The restriction of this function on the constraint surface is E l \c= 0. The action of the 
Hamiltonian vector field X^c on E 1 is 

X dC (E 1 ) =: E 2 = -ujxq sin (ut) + — cos (ut) - — , (49) 

m m 

and the restriction of E 2 on the constraint surface is E 2 \c = 0, and more important, the 
equation E 2 \c= is precisely an evolving constant 

p = —mcux sin (ut) + p cos (ut) . (50) 



Note also that the action of X^c on E 2 gives again E 1 , and the process ends. In other 
words, the evolving constant (|50|) was obtained from the application of the Hamiltonian 
vector field X^c on E 1 , and viceversa. 

The full solution. In the present case the constraint surface is coordinatized by the 
coordinates of the physical phase space (x ,p ) and by the internal time t. Therefore, 
Eqs. (0), and (0) acquire the form 

x = X(t; x ,p ) = x cos (ut) + — — sin (ut) , 

mu 

t = T(t; x ,p ) = t 

p = P(t; xo, po) = —muxo sin (ut) + po cos (ut) , 

Po 2 1 

p t = P T (t; x ,p ) = -mu 2 x 2 , (51) 

2m 2 



Of course, last equations are also (0), and (|20|). Notice that Eqs. (§), and (§) acquire 
the form 

x = cos (ut)x sin (ut)p , 

mu 

Po = cos (ut)p + mu sin (ut)x , (52) 
and the dependence of the orbit coordinate x 1 = t, see QlCf) , is 



t = T(x,t,p,p t ) =t. (53) 
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Last equations constitute the full solution of the classical dynamics of the system. 



Notice that the internal time variable x 1 = t = T(x,t,p,p t ) — t is not a physical 
observable because the Poisson bracket with the first class constraint does not vanish. 
Nevertheless, when we take the full solution into account we can express t = T(x,p,p , x ), 
given by 



(±muxx + ±ppo) 
cosut=- — -, (54) 



with H = ^Po 2 + \rnujXQ 2 . The above expression is an evolving constant. From this 
point of view, the internal clock t defines a two-parameter family of physical observables 
on the physical phase space; x, and p being the parameters. So, the internal clock t 
becomes a physical clock t(x,p), namely, a physical observable when the full solution is 
considered. We restrict the analysis to a branch of the above multivalued function to 
compute the time t(x,p) at which the particle reaches the position x and the momentum 
p evolving from an initial position x and momentum p Q 

,( \ v 1 (\muxxQ + ^pp \ ,„v 
t(x,p) = T(x,p; x ,po) = -arc cos A . (55) 

u \ H J 



Or in terms of x only 
t±(x) 



( 



-arc cos 



\mojxxQ ± 



2^ (^o 



Ha 



(56) 



These classical expressions have a quantum version as we will see later. 



B. Quantum dynamics 



At quantum level, as Dirac showed, the physical states are those killed by the first class 
constraint. We associate abstract operators with the classical coordinates and momenta, 
given by 

x^X , t^f , p^P , pt-^Pr, (57) 



which satisfy the Dirac rule 
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[X,P\=i% , [f,Fr]=ih, 



(58) 



and by inserting these operators in the quantum constraint C \ ip), this equation becomes 

( P T + 7^ + \^xA | tP) = . (59) 



Any physical state can be expressed in terms of the single quantum number of the har- 
monic oscillator, in abstract Dirac notation 

n 

f=J2\n){n\ . (60) 



In last expression, the physical states | ip) are 'frozen' (i.e. they are abstract vectors), the 
complex coefficients C n are constants. Notice that we have not choosen the coordinate 
basis yet. Taking a 'coordinate representation' | x,t) where the operators acquire the 
form 

(x,t \ X \ ip) = x{x, t | ip) , 
(x, 1 1 T | if)) — t(x, t | if)) , 

(x,t\P\if)) = j^(x,t\iP), 

(x,t\P T \if)) = ^(x,t\if)), (61) 

any physical state vector | if)) is expanded in the coordinate basis | x, t) as 

(x, t\if)) = ifj(x, t) = C n (x, t\n) , 

n 

(x',t f \x,t) = J2( X '^' \ n )( n \ > ( 62 ) 



with (x,t | n) = e~% Ent f n (x), E n = fiu)\n+^). Thus, in the Dirac framework, the 
coordinate representation is nothing but the 'Heisenberg picture' of the standard quantum 
mechanics, where the coordinate basis | x, t) is 'rotating' and the physical state | if)) is 
fixed (see Eq. flBUD where the coefficients C n are constant complex numbers). 



Schrddinger equation. In addition, we can build a 'Schrodinger basis' from the 'Heisen- 
berg basis' | x, t). In this 'Schrodinger basis', which we denote by | x), the state vector is 
'moving around' the 'fixed basis' | x). Explicitly, 

if>{x,t) = (x | ip(t)) , Schrddinger basis | x) , (63) 
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with 



\m)=T, C n(t)\n) = Y,C n e 

n n 

(x | n) = f n (x) . 



n) , 



(64) 



Taking the derivative with respect to the coordinate t of | ip(t)), the familiar Schrddinger 
equation emerges in the formalism [T^ 



ihj t i m) = h i m) 



(65) 



with H = ^P 2 + \muj 2 X 2 . As usual, the physical vector | if)(t)) evolves in t while the 
coordinate basis | x) is fixed. In other words, if we consider the system composed of the 
harmonic oscillator plus the clock together, we are describing the evolution of the degrees 
of freedom of the harmonic oscillator with respect to the internal clock itself, that is to 
say, the evolution of one part of the system with respect to the rest of it. In the next 
section, we will carry out the same procedure we applied here in order to analyze the 
meaning of evolution in generally covariant quantum theories. 

Quantum evolving constants. Let us now go to the quantum version of the evolving 
constants. The Hilbert space is built with the implementation of the physical state vectors 
I "0) = J2 n Cn I n ) in the reduced Hilbert space Ti r associated with the physical phase 
space of the harmonic oscillator. In the present case 



1p(x Q ) = (X | ijj) = ^2C n f n (x Q ) 



(66) 



The inner product in 7i r 



= / dfx(x ) i)*(x )(p(x ) 



(67) 



can be determined with the condition that the operators xq, and po be hermitian opera- 
tors. 



Thus, in the classical expression 



x = X(t; x ,p ) = xq cos (uot) + -^L sin (cut) , 

moj 



(6* 



x , and p are physical observables given by the Eq. ([54]) and they become operators 
acting on T~C r , so the quantum version of the classical evolving constant is 

x(t) = x(t; Xo,Po) = Xq cos (u>t) + — — sin (ut) , (69) 

rrujj 
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which is the well-known evolution equation for the position operator X in the Heisenberg 
picture. 

In addition, the classical expression 

p = P(t; x ,p ) = —muxQ sin (u>t) + p cos (uxt) , (70) 

has its quantum analog 

p(t) = p(t; xq,po) = —mojxo sin (ut) + p cos (out) , (71) 



and finally 



p 2 1 

p t =Pt(t;x ,p ) = --2- - -mu 2 xl. (72) 



In summary, for the case of parameterized systems, the quantum version of the evolving 
constants equations constitutes the Heisenberg equations for the physical operators in- 
volved in each particular theory. In the case of the harmonic oscillator, Eqs. (|69| ) and 

(□)• 



Time operator. The classical expression ( p5|) becomes an operator T(X, P) = 
t(X, P;xo,po) which is defined on the reduced Hilbert space 7i r . Taken arbitrarily the 
order of the operators, we have 

T(X, P) = -arc cos M L 2m — J , (73) 



with Hq = 2^;j?o + \m(jjx\. From this operator, we can compute the 'time of arrival' 



operator T(X) 



T±(X) = —arc cos 

UJ 



(74) 



associated with the time at which the harmonic oscillator is detected with an apparatus 
located in x = X. The 'time of arrival' operator for a free particle has been studied in 
T3fl . The analysis of the 'time of arrival' operator for the harmonic oscillator deserves to 



be studied. 
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III. SL(2,R) MODEL WITH TWO HAMILTONIAN CONSTRAINTS 



A. Classical dynamics 



Let us see how the relative evolution looks in a non familiar generally covariant model. 
A nonlinear generally covariant system with two Hamiltonian constraints and with one 
physical degree of freedom was introduced in [|TJJ]. This model mimics the constraint 
structure of general relativity. Here, we continue the study of this model. In particular, we 
display the full set of evolving constants required in its classical and quantum dynamics. 
Moreover, for a Schrodinger-like equation of motion arises in its quantum dynamics, we 
compare the meaning of time (evolution) in both, evolving constants and Schrodinger-like 
equation, viewpoints. 

First, a brief summary of its classical dynamics, for more details and its physical 



interpretation see Ref. [14]. The model is defined by the action 



S[u, v, N,M,\} = l[ dt[N(Vu 2 + v 2 ) + M (Vv 2 + u 2 )], (75) 

where 

Vu=^{il-\u), Vv=jj(i!+\v); (76) 

the two Lagrangian dynamical variables u = (it 1 , u 2 ) and v — (v l , v 2 ) are two-dimensional 

real vectors; N, M and A are Lagrange multipliers. The squares are taken in R 2 : u 2 = 
u ■ u = (u 1 ) 2 + {u 2 ) 2 . The action can be put in the Hamiltonian form 

S[u, v,p,n,X m ] = [ dr • p + €?• 7? - A m C m l . (77) 



The canonical pairs that coordinatize the unconstrained classical phase space are (u 1 ^ 1 ), 
(u 2 ,p 2 ), (v 1 ,^ 1 ), and (v 2 ,ir 2 ). Also A 1 = N, X 2 = M, and A 3 = A. The first class 
constraints have the form 

C 2 = l(i? 2 -u 2 ), 

C3 = u ■ p — V ■ 7T , (78) 



which algebra is isomorphic to the sl(2, R) Lie algebra 
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{Ci, C2] — C3 
{C 1 ,C 3 } = -2C 1 

{C 2 ,C 3 } = 2C 2 . (79) 



The classical dynamics is the unfolding of the gauge symmetry generated by the Hamil- 
tonian vector fields 

x dc\ = ~P- V u -v-V n , 
XdC 2 = -7T • V„ - u ■ Vp , 

-XdO, = -u- V u + v- V v + p- V p -Tf- V,, (80) 



associated with the first class constraints of the model. 

The physical phase space can be coordinated by the points (J, <fi, e, e'), and these 
physical observables have the following form 



1 2 

up - 


1 2 
- p u 




6 = 1 1 2 

\u p - 


- p 1 u 2 \ 


1 


,_ vrV- 












\ix l v 2 - 


- w 1 ^ 2 


1 


J = \u l p 2 - 


-pV| 




= arctan 


u x v 2 — 


P^TX 1 


u l v l — 


p l 7T l 





The Poisson brackets between J and <p in the reduced phase space reads 

{J,0} = ee'. (82) 



Classical evolving constants. Finally, the relation between the Lagrangian variables 
[u, v) and the physical states (J, 4>, e, e') 



u v 1 + ee'u 2 v 2 



cos 6 + 



-u 1 w 2 — ee'u 2 v 1 



sin 1 



J, 



(83) 



which leads to the notion of evolving constants of the system [|]|| . The evolving constants 
give the evolution of the Lagrangian variables of the system in a gauge invariant way, i.e., 
for any fixed physical state of the system (J, 0, e, e'), Eq. (|83"D gives the change of one of 
the four coordinates as a function of the other three coordinates, say 

1/ -e'x(zcoS(j)- ysin0) + eJ 

C7 (ar, y, 2; J, 0, e, e ) = ? — — . (84) 

e[y cos <p + z sin 0J 
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This relative evolution among the coordinates is gauge invariant. In addition, for any 
fixed x, y, z last equation gives a three-parameter family of physical observables, the 
three parameters are the three coordinates x,y,z, on the physical phase space. 

Generation of evolving constants. We start with the evolving constant (jS3|) , and define 
the evolving function E 1 



E l (u,v,p,Ti) :-- 



u v + ee'u 2 v 2 



cos 6 + 



u x v 2 — ee'v^v 1 



sin <b — J . 



15) 



The restriction of E 1 on the constraint surface vanishes, E 1 \c= 0. The action of the 
Hamiltonian vector field X d H 1 on E 1 is 



XdH^E 1 ) =: E 2 



Piv 1 + ee'p 2 v 2 



cos — 



Piv 2 — ee'p2V 1 



sin i 



(86) 



and the restriction of E 2 on the constraint surface vanishes, so E 2 \c= gives the evolving 
constant 



Piv 1 + ee'p 2 v 2 



cos 6 + 



Piv 2 — ee'p2V 1 



sin 6 = 0. 



*7) 



The action of X^Hi on E 2 gives zero, so the process ends. Now, we compute the action 
of the Hamiltonian vector field on E 1 



XdHiiE 1 ) ='■ E 3 



U 7Ti + ee'u 2 TT2 



COS <t> + 



u l H2 — ee'u 2 Tii 



sm i 



and the restriction of E 3 on the constraint surfaces vanishes, so E 3 \c= gives the 
evolving constant 



u l iti + ee'u 2 iT2 



cos 6 + 



v}ii2 — ee'u 2 7Ti 



sin 6 = 0. 



?9) 



The action of XdH 2 on E 3 gives zero, so the process ends. Finally, the computation of 
the action of the Hamiltonian vector field Xjd on E 1 



X dD (E l )=: E* = -E l -J, 



(90) 



so we recover the original evolving constant we start with, and no more evolving can be 
obtained from 



The full solution. Eqs. (|T9"|), and ( p0|) acquire the form 

u 1 = u\u 2 yy-, j,<j>,e,(?) 



—e'u 2 (v 2 cos 6 — v 2 sin 6) + eJ 



e(v l cos 6 + v 2 sin . 
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u 2 


= u 2 > 


(u 2 . 






J,q 




e' 


v 1 


= V 1 


{u 2 


v 1 


v 2 - 


■ I 6 




e ' 


v 2 


= v 2 


(u 2 


V 


,v 2 : 


; J,<f 




y 


Pi 


= P1( 


u 2 , 


v\ 


v 2 ; 


JA 






P'2 


= p 2 ( 


■ 1 

u 2 , 


1 

v\ 


2 


J A 




/\ 


7Tl 


= u v 


[u 2 ., 






J A 




e' 


7T2 


= n 2i 


[u 2 . 


v\ 


v 2 ; 


J A 


>,e, 


< 



V, 

v 1 , 
V 2 , 

e (v 1 sin cf) — v 2 cos c/>) , 

e (v 1 cos + v 2 sin 

ew 2 ?; 1 + e'Jsin 
(v 1 cos + v 2 sin </ 

ew 2 f 2 — e' J cos 
(w 1 cos + v 2 sin <7 



(91) 



and the Eqs. (|3|), and (0) are precisely the Eqs. flgip while the Eqs. QTD|) acquire the form 

u 2 = TPiu^v^pititi) = u 2 , 

v 2 = V 2 (u\v\ Pi ,7T l )=v 2 . (92) 

So, the dynamics of this model can be described in a relational fashion way. The differ- 
ence with respect to parameterized systems, as the example of the harmonic oscillator 
previously analyzed, is that in the present case a single internal time variable is not 
enough, rather, we need three internal time variables. In the way we have expressed the 
full solution (|9~l]), vt 2 ,v x ,v 2 are clocks, i.e., once the component u 2 of the position of the 
first particle, and the position (v 1 , v 2 ) of the second particle are known, the change of the 
component u 1 of the first particle and the change of the momenta of both particles p, n 
are also known when the system is an particular physical state (J, 0, e, e'). Therefore, the 
full relational evolution of the system is expressed in terms of three internal clocks u 2 , v 1 , 
and v 2 . 



B. Quantum dynamics 

At quantum level, the model is characterized by the following set of observables 



J 


to, 


e, 


e'> 


= mh to. 




e 


to, 


e, 


e') 


= e to, e, 






TO, 


e, 


e') 


= e to, e. 


y> 



and the physical states are given by 

IV>>= £ C m , e ^ |m,e,e'), (94) 

me,e' 
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in abstract Dirac notation. In the 'coordinate representation' | u, v, a, (3), which is nothing 
but the Heisenberg picture in standard quantum mechanics because all the coordinates 
(u, v, a, (3) are put at the same level, the state reads 

i()(u, v, a, (3) = (u, v, a, (3 \ ip) = Cm,e,e'{u, v,a,(3\ m, e, e') , (95) 

mf,(' 



with (u,v,a,f3 \ m, e, e') = e im ^ a e//3 )j m (^). Thus the basis | u,v,a,f3) is 'rotating' 
and the state | if)) is fixed, i.e., the coefficients C mee i are constant complex numbers. The 
'coordinate representation' appears as the most 'democratic' basis because it does not 
prefer one coordinate more than the others. 

Schrddinger equation. In the same sense that in parameterized systems we were able 
to build a 'Schrddinger basis' from the Heisenberg basis, we can do the same here, and 
rewrite the physical state (0). In the present example, we can build two Schrddinger 



bases | u,v,@), and | u,v,a). In the first one, the physical state vector (|94|) is expressed 
as 

ip(u, v, a, /3) — (u, v, j3 | i(j(o()) , (96) 

with 

| </>(<*)> = ]T C m ^{a)\ m^e') = £ C m ^e^ a \ m^,e>) , 

me,e' me,e' 

(u, v, 0\ rr^7, e>) = e~ m ^J m . (97) 



Taking the derivative with respect to the coordinate a of | i/j(a)), a Schrddinger equation 
emerges in the formalism 

i h * | 4(a)) = _£o 3 4 | V(a)) , (98) 
da e 

and the physical observable O34 has the form 

(u, v, (3 I 34 I V) = (99) 

in the 'Schrddinger basis' | u,v,(3). As expected, in the Schrddinger basis | u,v,/3) , the 
state I "0(a)) evolves while the basis | u,v,(3) is fixed with respect to a. This is not 
a matter of terminology, in fact the evolution equation (^) is well defined, and we are 
really able of describing evolution under this picture, namely, to describe the change of 
the some part of the whole state with respect to rest of it, in complete agreement with 
the spirit of relationism. 
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In the second Schrodinger basis | u,v,a), the physical state vector ([54]) is expressed 

as 

t/j(u,v,a,/3) = (u,v,a\ip(0)), (100) 

with 

me,e' me,e' 

(u, v, a\ mT^ e') = e im£Q J m . (101) 

Taking the derivative with respect to the coordinate (3 of | 4>(j3)), a Schrodinger equation 
emerges in the formalism 

in lp 1 m) = l dl2 1 m) ' (102) 
and the physical observable Ou has the form 

i da 



% d 

(u, v, a | O12 | i>) — - — (u, v , a \ tp) , (103) 



in the 'Schrodinger basis' | u, v,a). 

Quantum evolving constants. The quantum version of the evolving constants is as 
follows. More precisely, the quantum version of the full classical solution (|91~D is expressed 
as 

-e'eiv 2 cos 6 — t/sin 6) + J 
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= u l 


in 2 


t/ 1 
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sin q 


), cos 4 
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>, cos 
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P2 


= P2 


y, 
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>, cos 4 
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7Tl 


= 7Tl 
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>,e 
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7?2 


= 7T 2 


(u 2 


t, 1 


v 2 ; 
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sin </ 


>, COS (j. 


»,£ 


?) = 



v 1 cos — v 2 sin (p 
t (V sin cf> — v 2 cos cj)j 

e (Vcos0 + v 2 sm(pj , 

ev^v 1 + e 1 Jsin 
v 1 cos + w 2 sin ' 
eu 2 v 2 — e'Jcoscfi 
v 1 cos d> + i> 2 sin 6 



(104) 



The meaning of the first equation in ( |104| ) is the following: we have to take the mean 



value of the operator u 1 with respect to generic states | ip) of the reduced Hilbert space 



H r of the model |U| 
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In summary, the quantum dynamics of parameterized systems can be described in 
terms of a 'Schrodinger equation' or in terms of the 'Heisenberg picture'. The Schrodinger 
equation arises as a consequence of the Dirac quantization, as we have seen for the case of 
the harmonic oscillator. On the other hand, in the SL(2, R) model, we were able to build 
two (dependent) Schrodinger equations, and thus to identify two (dependent) internal 
time variables a and (3 with respect to which the physical states of the SL(2, R) model 
evolve. This does not mean that is always possible to single out in general an internal 
time variable, given by a Schrodinger equation, in generally covariant theories once the 
Dirac quantization has been performed. Therefore, in general, a Schrodinger equation 
does not arise in the formalism. The Schrodinger picture, when this picture emerges 
in the formalism as a consequence of the Dirac quantization, singles out one internal 
clock only. More important, the quantization of generally covariant theories based on the 
reduced Hilbert space (generalized Heisenberg picture) need M internal clocks, where M 
is the number of first class constraints. In the case of the SL(2, R) model, the clocks 
are u 2 , v 1 , and v 2 in the generalized Heisenberg picture. In the Schrodinger picture, the 
internal clock is given by a (or (3). 

Classical limit. Now, we compare the quantum evolving constants of the SL(2, R) 
model with those of the harmonic oscillator in order to get insights on the classical limit 
of generally covariant theories, and in particular of the SL(2, R) model. We expect that 
the classical limit of generally covariant theories should be attached to the concept of 
coherent states as it happens in standard quantum mechanics (parameterized systems). 
In the case of the harmonic oscillator, the coherent states are roughly those states | ip) in 
the reduced Hilbert space TC r such that the mean values (if) \ x(t) \ if)), and (if) \ p(t) \ ip) 
reproduce the classical behavior of the system. Of course this condition is not enough to 
single out the coherent states of the system. In addition, those states have also to minimize 
the uncertainty relations of position and momentum. Of course, these two conditions are 
still not enough to identify the coherent states due to the fact that both conditions are 
satisfied by both squeezed and coherent states. In the case of the parameterized harmonic 
oscillator a mechanism that identifies the coherent states is available following standard 
methods. It is natural to expect that a combination of the coherent states approach to 
the quantization of generally covariant theories [|TjJ with the full set of evolving constants 



of motion required in their quantum dynamics displayed here could bring the classical 
limit of constrained systems. 



IV. CONCLUDING REMARKS 



We have displayed the full solution of the relational evolution of the degrees of freedom 
of fully constrained theories with a finite number of degrees of freedom (see Eqs. (|T5|), 
and (0)). Our procedure follows from the embedding equations of the coordinates and 
momenta in the unconstrained phase space (see Eqs. (D, and ©) plus the expressions 
of the M internal time variables (see Eq. (|H])). The form of the solution containts 
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all the evolving constants of motion needed in the description of the classical dynamics 
of fully constrained theories, i.e., we have given the full mathematical solution to the 
Rovelli's point of view on the 'problem of time' pioneered in Refs. 0-§|. Of course, the 
physical (and phylosophical) interpretation is due to Rovelli. Also, we have explored 
a method to generate those evolving constants. This method consists in the repeated 
application of the Hamiltonian vector fields associated with the first class constraints on 
some initial evolving constant. Combining the expressions of this evolving constants with 
the expressions of the physical observables the full relational evolution of the coordinates 
and momenta is obtained. Finally, we have also analysed on a general setting the quantum 
version of the relational evolution of the degrees of freedom of fully constrained theories. 

To find the full solution of the relational evolution of the degrees of freedom for gravity, 
matter fields coupled to gravity (see |TB[ for the first steps), topological quantum field 



theories, or for a background-independent string theory constitutes one of the challenges 
of the new millenium. 
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